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ABSTRACT 

Residual-gas expulsion after cluster formation has recently been shown to leave an 
imprint in the low-mass present-day stellar mass function (PDMF) which allowed 
the estimation of birth conditions of some Galactic globular clusters (GCs) such as 
mass, radius and star formation efficiency. We show that in order to explain their 
characteristics (masses, radii, metallicity, PDMF) their stellar initial mass function 
(IMF) must have been top-heavy. It is found that the IMF is required to become 
more top-heavy the lower the cluster metallicity and the larger the pre-GC cloud- 
core density are. The deduced trends are in qualitative agreement with theoretical 
expectation. The results are consistent with estimates of the shape of the high-mass 
end of the IMF in the Arches cluster, Westerlund 1, R136 and NGC 3603, as well 
as with the IMF independently constrained for ultra-compact dwarf galaxies (UCDs) . 
The latter suggests that GCs and UCDs might have formed along the same channel or 
that UCDs formed via mergers of GCs. A fundamental plane is found which describes 
the variation of the IMF with density and metallicity of the pre-GC cloud-cores. The 
implications for the evolution of galaxies and chemical enrichment over cosmological 
times are expected to be major. 

Key words: stars: formation - stars: mass-function - stars: early-type - stars: late- 
type - globular clusters: general 



1 INTRODUCTION 

1.1 Are PDMF variations evidence for IMF 
variations? 

The stellar present-day mass function (PDMF) is ob- 
served to be different for individual globular clusters 
(GCs) of the Galaxy. Their low-mass stellar content differs 
siRnificantly (iMcClure et al.l Il986l. iDiorgovski et all 1 19931 : 



iDe Marchi et al.ll2007l : |Paust et alll2010l ) 



These differences are claimed to be explainable by 
a universal initial mass function (IMF ) and secular two - 
body relaxation driven evolution alone (|Leigh et al.ll2012h . 
in agreement with earlier work that shows that an ini- 
tially rising IMF with decreasing stellar mass ca n bend over 
in the course of a Hubble t i me llVesperini fc H cggic I997I; 
iBaumgardt fc Makind 120031 '). iDe Marchi et al.l (|2007l ') how- 
ever noticed that clusters which are strongly depleted in 
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low-mass stars have a low concentration while no such clus- 
ter with a high concentration exists. This trend has been 
ar gued to be incompatible with standard secular evolution 
by iDe Marchi et all ll200'itl unless the majority of clusters 
with a flat or eve n declining low-ma ss PDMF are post-core 
collapse clusters. IPaust et al.l ()201Clt l question the existence 
of such a trend since no very strongly depleted cluster is 
present in their sample. Their data however sti l l close ly fol- 
lows the relation proposed bv lDe Marchi et al.l |20o3). 

Primordial mass-segregation appears to come to the res- 
cue of the secular evolution picture, since if low-mass stars 
form preferentially at the cluster outskirts they are easily 
stripped off the cluster by the tidal-field of the host galaxy. 
There is then no need for the cluster to go into core-collapse 
which would otherwise be necessa ry to drive the preferen - 
tial evaporation of low- mass stars l|Baumgardt et al]|2008l ). 
However, the se authors al s o sho w that some clusters re- 
ported by [De Marchi et al.l (j2007l ) which are most strongly 
depleted in low-mass stars, can only be understood if clus- 
ters fill their tidal-radii. Tidally underfilling clusters can not 
reach the observed degree of low- mass star depletion, eve n 
when they are close to dissolution (|Baumgardt et al]|2008l ). 
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Young star clus ters of the Galaxy a re however observed to be 
very compact (|Lada fc Ladalfioosi ). not closely filling their 
tidal limits and such initial conditions cannot be the domi- 
nant cause of low-mass star depletion. 

Is the concentr ation-PDMF trend observed by 
iDe Marchi et al.1 (|2007l l therefore the first-time evidence for 
a variable low-mass IMF? iMarks fc Kroupal (|2010h show 
that clusters wit h a flat or declining PDM F have the largest 
metallicity in the lPe Marchi et al.l (j2007l ) sample of clusters. 
Such a trend is difficult to understand in the dynamical evo- 
lution picture since it is unclear how dynamics could possi- 
bly know about the metal-content of the cluster. The trend 
is also difficult to understand if the low-mass IMF were to 
vary since relatively fewer low-mass stars would have to form 
with decreas ing metallicity, in contradiction to theo retical 
expectation jBastian et all boid : iKroupa et al] 1201^ 1 . Fur- 
thermore the concentration-PDMF trend is also difficult to 
reconcile with the assumption of a varying low-mass IMF: 
In dense, i.e. highly concentrated clusters, low-mass stars 
might disappear through merging in dense clusters so that 
low-mass stars would be expected to be underabundant in 
highly concentra t ed clu sters, which is again contrary to the 
iDe Marchi et al.1 (|2007l ) observation. 

Thus, if standard two-body relaxation alone is not suf- 
ficient in removing low-mass stars and a varying low-mass 
IMF appears not to be a feasible solution, whe re do the ob- 
served differences in the PDMF come from? iMarks et al.l 
l|2008t ) propose a residual-gas expulsion scenario which 
adresses all these issues at once. They show that quick gas 
removal from compact, primordially mass-segregated clus- 
ters with binaries leads to low-concentration clusters which 
are depleted in low-mass stars, while initially binary-rich 
clusters w ith slow gas remo v al ret ain their input IMF, as ob- 
served bv lDe Marchi et al.1 |20o3). The expansion following 
gas expulsion (accompanied by the loss of primarily low- 
mass stars) naturally leads to tidally-fiUing clusters with 
low concentration. This process is expected to be metal- 
licity dependent since metal-rich material will couple bet- 
ter to the rad iation driving the removal of the residual-gas 
l|Marks fc Kro upa 2010). The resulting quicker gas expul- 
sion, in turn, leads to more metal-rich clusters being more 
strongly depleted of low-mass stars, as observed. This gas 
expulsi on scenario is i n qua l itativ e agreement with the find- 
ing by IStrader et al.l (|2009l . I2OIII ) that mass-to-light ratios 
for M31 GCs are lower in relatively more metal-rich envi- 
ronments. 

Thus, invoking gas expulsion there is no need for a varia- 
tion in the low-mass IMF. In contrast, for stars more massive 
than « IMq theoretical arguments suggest that the massive 
star content should depend on the ambient conditions. 

The Jeans-mass (Mj oc p~^/^r^/^) is higher for de nser 
l|Larsonlll998l :lBate fc Bonncll"2005l: lBonnelI et al.l2006l ') and 
warmer gas (Klcsscn ct al. 200^ such that a higher av- 
erage mass for stars is expected when favourable condi- 
tions are met. As both the temperature and density of the 
gas depend on its metallicity through less efficient cool- 
ing in lower metallicity environments, more massive stars 
should form from low metallicity gas. Additionally, form- 
ing stars self-regulate th eir masses via radiative feedback 
(jAdams fc Fatuzzd[l996l ). Feedback is also expected to be 
metallicity dependent, as photons couple less efficiently to 
gas of lower metallicity. Anew, a lower metallicity should 



favour the formation of more-massive stars. Therefore both 
the Jeans-mass- as well as the self-regulation-arguments 
suggest the fraction of high to low mass stars to increase 
with decreasing metallicity. 

Also, star clusters have been found t o be rather com- 
pact when they reach their denses t state (iTesti et al.l 19981; 
Kroupal l2005l : IScheepmaker et aL I I2OO9I : iMarks fc Kroupal 



The resulting density of stars is expected to have an 
influence on the IMF (|Bonnell et a"l][l998, : ,Shadmehri.2004l ') . 
lElmegreen fc Shadmehril (|2003l ') concluded that densely 
packed stars should lead to top-heavy IMFs, i.e. an IMF 
wit h more massiv e stars than expected from the canoni- 
cal iKroupal l|200ll ) IMF (equation [2] below), in the most 
massive cluste r s (e.g . GC progenitors) through merging. 
IWeidner" et al.l {20m show, using a geometrical argument, 
that a high density will be a problem for the formation of in- 
dividual stars in clusters more massive than ~ 10^ Mq. And, 
cluster formation at high star formation rates may lead to 
top-heavy IMFs by the heating of molecular gas through 
superno vae generated cosmi c rays (the Papadopoulos CR- 
heating, jPapadopoulosI I2OI0I ') . 

Finding observational evidence for a varying high-mass 
IMF is difficult. In GCs, stars more massive than IMq have 
long since evolved away from the main-sequence and cannot 
be observed. Rather weak direct observational evidence for 
top-heavy IMFs in young starburst clusters exists: 

The Arches-cluster (2 — 4 Myr), only ~ 25 pc in pro- 
jected distance from the Galactic centre, should be sub- 
ject to a very strong tidal-field and is a prime candidate 
which has been reported to host a non-canonical popula- 
tion of high-mass stars. The measureme nts however differ 
widely (03 = 1.6 5 — 2.3 in equati o n [2] Figer et al. 19991 : 



20091 ). 



IStohe et al.1 12OO2I : iKim et al.1 12OO6I : lEspinoza et alT 
The same holds true for the NGC 3603 cluster (2-3 Myr), 
about 8 kpc from the Galactic centre and 6 kpc from 
the Sun. Its reported MF slope ranges from 03 = 1.9 to 



2.3 (iNiirnberger fc Petr-Gotzensll200a:ISung fc Bessellll2004l : 



IStolte et all I2OO6I : iHaravama et al.1 l2008h . The most mas- 
sive known young cluster in the Galaxy Wester lund 1 (Wdl, 
4—7 Myr) in contrast h as as — 2.3 between .75 and 1.5 pc 
and 03 = 1-6 closer in jBrandner et al.ll2008l '). The 30 Dor 
region (~ 3 Myr) in the Large Magellanic Cloud (LMC), 
having a mass comparable to Wdl, also exhibits no de via- 
tion from the canonical IMF (jSelman fc Melnickllioosl . ex- 
cludi ng the central cluster R136). R136 has Q3 = 2.2 — 2.6 
( Br andl et al.lll996l : iMassev fc Hunteijll998l : lAndersen et~aLl 
|2009j). All reported clusters however show evidence for mass 
segregation, i.e. the r esults might be biased throu gh dynam- 
ical effects (see, e.g.. jPortegies Zwart et aLll2007^ . 

Better evidence for top-heavy IMFs is seen in ultra- 
compact dwarf galaxies (UCDs). These are objects simi- 
lar to GCs, given the ideas of their formation from star 
cluster complexes ( Fellhaucr fc Kroupa 2005) or as the 
most massive GCs ((Mieske et al.l I2OO2I ) and are thus ide- 
ally suited to be compared to GCs. Observed dynamical- 
mass-to-light (Mdyn/i)-ratios appear too large when they 
are compared to expectations from canonical stellar popula- 
tion m odels with the canoni c al IM F (Dabringhaus en et al.l 
I2OO8I ). iDabringhausen et al.1 (I2OI0I ) explain these observa- 
tions by invoking top-heavy IMFs for them, i.e. UCDs had 
non-canonical IMFs above IMq. The subsequent, primar- 
ily gas expulsion driven evolution then turns them into ob- 
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jects wh ich resemble the obse r ved p roperties of UCDs to- 
day. The lDabringhausen et all (|2010l ) simulations show that 
top-heavy IMFs with 1 < as < 2.3 for proto-UCDs are in- 
deed viable to explain the presently observed Afdyn/i'-ratios 
in them and their models resemble present-day UCDs in 
mass and size. This result is significantly strengthened by 
the overabundance of X-ray-bright UCDs which is in excel- 
lent agreement with the surplus of ne utron-stars and black 
holes expected for a top-heavy IMF (jDabringhausen et al.l 
I2OI2I ). 

In this contribution it is shown that evidence for a 
top-heavy IMF also emerges for GCs if residual-gas expul- 
sion from mass-segregated clusters starting with a universal 
low-mass IMF is r esponsible for the obser ved low-mass star 
depleted PDMFs jPe Marchi et"al]|2007f ). as suggested by 
iMarks et al.l l|2008l l . In order to deplete a universal low-mass 
IMF in GCs to the observed degree through the gas throw- 
out phase the IMF must have been top-heavy in order for 
the high mass stars to be sufficiently destructive to the gas 
content of the forming GCs. 



1.2 Tracing the high-mass IMF in GCs via 
residual-gas expulsion 

The IMF at the high-mass end can never be directly seen 
in Galactic GCs, since massive stars are short-lived while 
GCs are old objects. Stars above the turn-off mass for GCs 
(« O.8M0) cannot be observed. However, the distribution of 
high-mass stars, which were present in a GC after its birth, 
can possibly be traced by indirect means. 

Young gas-embedded star clusters expell their left-over 
gas from s tar formation due to stellar feedback over a time- 
scale tm (|Kroupa et al.ll200ll ). For a given mass, size and 
star formation efficiency of a GCs pre-gas expulsion cloud- 
core (containing stars-fgas), tm is determined by the num- 
ber of massive O- and B-stars, which are the main donators 
of energy into the gas and, in turn, drive the gas expulsion 
process. Thus, the shape of the high-mass end of the IMF 
of stars defines tm- For a star formation efficiency (SFE), 



(1) 



0.1 < e < 0.5, rapid gas expulsion (short tm) gives 
rise to a phase of strong expansion following gas expul- 
sion and the cluster will loose sta rs over the tidal bound- 
ary l|Baumgardt fc Kroupal |2007| ). If a clus ter is mass- 
segre g ated at birth, as i ndica.t e d by theory (iMarks et al.l 
I2OO8I : iBaumgardt et all I2OO8I: lAllison et al.l |2009|) and 



observations |Littlefair et alj 120031 : iMcCradv et al.l I2OO5I : 
IChen et al ]|2007l ). the cluster will preferentially loose its low- 
mass stars. The form of the high-mass IMF can thus leave 
an imprint on the low-mass PDMF via the gas expulsion 
process. 

Using the iDe Marchi et all l|200'iD concentration, c — 
\og{rt/rh), versus low-mass PDMF slope diagram as a di- 
agnostic tool, SFEs, e, pre-GC cloud-core masses, Md, half- 
mass radii, r^, and densities, pci, at star cluster birth are 
constr ained for the sample of 20 GCs in iDe Marchi et al] 
l|2007t l for which PDMFs have been measured down to 
0.3 M(^ by comparison of o bservations with A'^-body mod- 
els (|Marks fc Kroupal I2OI0I '). These cluster parameters will 



in the following be referred to as the initial conditions of 
GCs and are reported in Table [T] 

In order to remove the gas from these clusters com- 
pletely it needs to travel between ~ 0.5 — 1.5 pc from the 
cluster centre to leave the cluster (compare to the rh re- 
ported in Table [!}. If the gas travels with the speed of sound 
in the ionised ISM (about 10 km s"^ « 10 pc Myr"^) the 
gas expulsion time, tm, becomes 0.5 — 1.5 x 10^ yr. This is 
an estimate for the shortest timescale on which gas can be 
expelled from a star cluster. It will be shown that gas cannot 
be removed within tm ~ 10^ yr (= tm for a cluster of size 
1 pc) if the IMF above IMq has been canonical but that 
instead a top-heavy IMF in dependence of the environment 
needs to be invoked (as in equation [2] varies). 



1.3 The models 



The following analysis relies on the lMarks fc Kroupal (|201G| ) 
results. In order to obtain the initial conditions they as- 
sumed residual-gas expulsion to be the dominant phase 
driving the l ow-mass star dep l etion in GCs and to be the 
source of the lDe Marchi et al.l (|2007l ') concentration-PDMF 
slope relation. They used the residual-g a s exp ulsion A^-body 
model grid by IBaumgardt fc Kroupal l|2007h to infer the 
mass loss an d expansion through gas removal from embed- 
ded clusters. iMarks fc Kroupal 1 2010l) assume that post-gas 
expulsion cluster expansion is only weak compared to expan- 
sion driven by gas expulsion. Baumgardt ot al. (2008) calcu- 
lated initial masses for the [Pe Marchi et al.l (|2007f ) sample 
of GCs assuming two-body r elaxation driven low-mas s star 
depletion without gas, which iMarks fc Kroupal (|201G| ) used 
as cluster masses after residual-gas throw out in order to 
account for mass loss through long-term secular evolution. 
This finally allowed to trace-back present-day cluster masses 

and sizes to th eir initial values. 

While the iMarks fc Kroupal l|2010l ) results depend on 
the validity of these assumptions, it is the first time that 
initial conditions have been constrained for a large sam- 
ple of real clusters. Undoubtedly improved initial conditions 
could be obtained by running A-body models including gas 
expulsion and two-body relaxation driven evolution after- 
wards to self-consistently compare observational data with 
computations. Such A-body models of initially compact 
Mod > lO^A^Q clusters are currently however prohibitively 
expensive computationally. For the time being their results 
constitute the best available homogeneously obtainted con- 
straints on GC initial parameters. The results derived here 
should therefore not be considered to be a very final word 
but they allow at the very least an important insight into 
star formation in massive star-bursting clusters at high red- 
shift. 



2 CONSTRAINING as 

Observationally inferred stellar initial mass functions 
(IMFs) of stars are found to be indistinguishable from the 
canonical form, ^c{m). Stellar IMFs are conveniently repre- 
sented by a multi-part power-law that describes the number 
of stars, dN, forming in the mass interval [m,m -f dm], 



dN /dm — ^(m) = k Oi x m 



1,2,3) 



(2) 
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where fc is a normalization constant and the Ui warrant con- 
tinuity at the edges of the power-law segments. For main- 
sequence stars the canonical IMF has the slopes qi,c = 1.3 
for tu/Mq G [0.08, 0.5] and a2,c = = as.c = 2.3 fo r m/M@ G 
[0.5, 1] and [l,mmax], respectively (iKroupalbOOll ). The phys- 
ical stellar upper-mass limit is m max 5, ISOMp, whic h de- 
pends on the mass of the cluster (jWeidner et ahllioiol ). 

Assuming the IMF to be canonical below IMq , we seek 
to infer the slope, as, in order to remove the residual-gas 
within a given time-scale, tm. 

The idea is to compare the energy required to remove 
the gas, -Egg, with the energy, -Eqej' provided by O- and 
B-type stars within tm- Eq^ is the difference between the 
initial binding energy of the cluster (with gas) and the bind- 
ing energy after the gas has been expelled. 



^OB 



(3) 



For a mass, Mpi, distributed according to a Plummer density 
profile, the potential energy is 



^OB 



32 



r-pi 



(4) 



where rh = 1 .305 rpi is the characteristic Plummer radius 
(|Kroupall200l ). Then equation (O becomes 



E'^^ = 1.305 X 



1.305 X 



37r G I'M, 

32 
37r G 



i-h,f 



32 



Mecl Mil 



(5) 



where Med ~ eAfd is the mass in stars, rh.j = ''fa is the 
initial half-mass radius (both from lMarks fc KroupallioiOl ) 
and rhj is the final half-mass radius of the cluster after the 
gas has been expelled. The last equality holds if the gas is 
removed adiabatically, i. e. slow with respect to the clusters 
crossing-time since then (|Hillj|l980l ) 

r.J = r.,. — . (6) 
The crossing-times, 

for the investigated sam ple of GCs as calculated f rom the ini- 
tial conditions derived in iMarks fc Kroupal (|2010l ) are indeed 
much shorter than the here assumed residual-gas expulsion 
times. 

The rate by which radiative plus mechanical energy 
from all stars is deposited in the interstellar medium (ISM) 

is, 



E 



E,{m) £,{m) dm 



(8) 



O.OSAIq 



where £^{rn)drn is the number of stars in the interval [m,m + 
dm] and Et{m) is the total energy output by a single 
star of mass m. This energy can be calculated from (see 
iBaumgardt et al.ll2008l ) 

logio i^./ergMyr-^ = 50 + 1.72 (logiom/M0- 1.55) , (9) 

i.e. the contribution through low-mass stars is negligible. 
Within a gas-removal time-scale, tm, an amount of energy 
equivalent to E^^{a3) — E{ai) ■ tm is released by mas- 
sive stars. The provided energy depends on the shape of 
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Figure 1. Energy provided by O- and B-type stars within 10^ yr 
as a function of the high- mass IMF index, «3, for the Galactic 
GC NGC 6752 (solid curve). The horizontal dashed line indicates 
the energy, £^og, required to remove NGC 6752's initial gas mass 
(equation [Sj . For an IMF index 03 = 1.27, the energy input from 
massive stars within 10^ yr is sufficient to overcome i?oB' 



the stellar IMF (equation |8} and therefore on the choice of 
Q3 and mmax. The latter has been found to have a negli- 
gible influence on the results and is arbitrarily chosen as 
niniax = 12OM0. In order to evaluate equation ([8]), the co- 
efficients k and ai in equation ([2} are found by normalizing 
the canonical IMF to Med, 



Moel = 



m (,c{m) dm 



(10) 



By then changing Q3 individually for the initial condition 
of each cluster (as constrained by iMarks fc Kroupal I2OI0I ) 
and calculating equation (|8]), a solution for as is obtained 
for which -Eqe ~ ^obi such that the energy is sufficient 
to remove the residual-gas. Using the same k and ai for all 
Q3 ensures that the number of low-mass stars, which provide 
the light after the massive stars extinguish, does not change. 



3 RESULTS 

Fig. [1] shows how the energy-input through massive stars 
changes with the IMF index as using the example of 
NGC 6752. The steeper the IMF (the lower 03), the more 
massive stars are available and the more energy is deposited 
in the ISM (solid curve) . NGC 6752 had to have a top-heavy 
IMF with as = 1.27 in order for the OB-star energy input 
to equal the energy required to remove the gas (dashed line) 
from its progenitor (-Bqb ~ ^ob)- 

The results for the estimated ci3-values in the sample 
of GCs are summarized in Tab. [1] For the respective initial 
conditions (Md, rn and e, also listed in Tab. [l| the results 
imply that the IMFs above IMq have been mostly top-heavy 
{as ~ 1 . . . 2.3, assuming the residual-gas is removed by ra- 
diation and winds from stars only). Slightly top-light IMFs 
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Table 1. The birth conditions for the GCs for which the PDMF has been measured from 0.8 down to 0.3 M(7) by[Pe Marchi et al. 
Column 1 and 2 list the cluster name and its metallicity. The next five columns give the constraints llMarks fc Kroupa 2010) on the 
inital GC cloud-core mass (Md), half-mass radius (r^,), density within the half-mass radius (pci) and lower limit to the star formation 
efficiency (eiow)- The last five columns are the assumed residual-gas expulsion time-scale, the here calculated required high- mass IMF 
slope {as) E'nd the corresponding number fraction (r/jvi equation 11611 and mass fraction (ryjv/i equation 1 1711 of stars more massive than 
IMq, respectively, and the total-mass fraction lost due to stellar evolution (/scv) over one Hubble-time. 
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6712 


-1.01 


2.17 


1.18 


0.16 


0.50 





1 


2.20 


0.10 


0.60 


0.50 


6752 


-1.56 


3.68 


0.24 


31.78 


0.15 





1 


1.27 


0.27 


0.96 


0.85 


6809 


-1.81 


1.99 


0.85 


0.39 


0.25 





1 


1.89 


0.13 


0.76 


0.65 


6838 


-0.73 


0.30 


0.47 


0.34 


0.50 





1 


2.60 


0.08 


0.44 


0.35 


7078 


-2.16 


17.30 


0.20 


258.13 


0.10 





1 


0.76 


0.56 


0.99 


0.89 


7099 


-2.12 


3.30 


0.23 


32.38 


0.15 





1 


1.29 


0.27 


0.95 


0.85 



are found for NGC 288, 6352 and 6838. A solution for Pal 5 
could not be found unless an extremely top-light MF is al- 
lowed (i.e. a canonical IMF more than suffices to remove the 
gas in Pal 5). This is likely connect ed to Pal 5 b eing close to 
dissolution, as already discussed in ; Marks fc Kr oupa (2010). 

To put these results into context, they are compared to 
independent theoretical and observational evidence for the 
existence of top-heavy IMFs in UCDs. Also, a comparison 
for different IMF observations in present-day Local Group 
starburst clusters (see Sec. [TJ is performed. 

For the starburst clusters their respective birth cloud- 
core masses and densities are needed for this comparison. 
These are here estimated using their present-day masses, 
half-mass radii and densities as compiled in Table [2] and 
assuming that the clusters have formed with e — 1/3. Given 
the youth of these clusters, no significant mass loss apart 
from gas-blow out is assumed to have occurred. Their cloud- 
core densities at birth, pci, were calculated similarly and 
assuming that their sizes did not yet change strongly during 
gas expulsion. The initial condi tions fo r UCD s are readily 
available from the Dabringhau sen et al.l l|2010l ) simulations. 

Fig.[2]shows the dependence of the high-mass IMF slope 
as a function of the pre-gas expulsion cloud-core mass. While 
the a3(Mci)-trend for UCD models and the constraints for 
GCs-|-starburst clusters are similar, both kinds of object ap- 
pear to form separate sequences running parallel to each 
other. However, when describing 0:3 as a function of the birth 
cloud-core mass-density, pci, i.e. incorporating information 
of the objects' sizes, the independently obtained data agree 
excellently with each other (Fig. [3]). 

Figs, m and [21 thus suggest that the high-mass IMF was 
more top-heavy (flatter) in more massive and denser envi- 
ronments. The expectation would indeed be that denser sys- 



Table 3. Coefficients (equation IllH for least-squares fits to the 
GC data (Figs. (2] [3] and |4]l. Additionally, coefficients for the de- 
pendence of Q3 on the stellar mass (M^ci = e X M^i) forming from 
the cloud-cores are given. Med Eind M^i are in units of 10® Mq 
and Pel is in units of W^Mq pc" ^. 



A 


P\ 


Ix 


■^lim 


Fig. 


logio A^ci 


-0.94 


2.14 


> 6.8 X IO^'Mq 


m 


logio A4cl 


-0.77 


1.59 


> 2.7 X IO^Mq 




logio Pel 


-0.43 


1.86 


> 9.5 X IO-'Mq pc-3 


m 


[Fe/H] 


0.66 


2.63 


< -0.5 


a 



tems form more massive stars with respect to the canonical 
IMF, if a channel for massive-star formation is the coag- 
ulation of proto-stellar cores since the collision probability 
is higher in denser sys tems, where the star s may also have 
lar ger accretion rates (|Murrav fc Linlll996l . Sec. [1} and/or 
the PapadopoulosI (|2010l ) cosmic ray heating is active in star 
bursts. 

Finally, Fig. [4] depicts 03 as a function of the present- 
day global metallicity for the GC and Local Group starburst 
cluster data. Metal-poorer environments appear to form flat- 
ter IMFs. Again, this reflects the theoretical expectation: In 
systems of low metallicity the Jeans mass is larger favouring 
the formation of m ore high-mass stars (|Adams &: Fatuzzol 
1 19961 : iLarsonl Il998l ) . A straight-forward least-squares fit to 
the GC data suggests the IMF to become canonical above 
[Fe/H] = -0.5. 

Thus, evidence for a systematic variation of the high- 
mass IMF m dependence of the birth environment emerges. 
In the following, each dependence of the suggested top- 
heaviness of Galactic GCs on their progenitor cloud-core 
properties is described by linear relations. Writing 



6 M. Marks, P. Kroupa, J. Dabringhausen and M. S. Pawlowski 



Table 2. Properties of starburst clusters in the MW and the LMC. From left to right the columns denote the cluster name, its age, 
stellar mass, Med, half-mass radius, r^, stellar density, Pecli within r^, and range of published high-mass PDMF slopes, 03, from the 
references given in the last column. 



Name 


age/Myr 


M,,i/Mq 


rh/pc 


Pcci/A^0 pc-=* 




References 


Arches 


2-3 


1.3 X 10* 


0.24 


1.12 X 10^ 


1.65-2.3 


1-6 


NGC 3606 


1-3 


0.7 X 10* 


0.2 


1.04 X IQS 


1.9-2.3 


3,4,7-9 


Wdl 


3-4 


5 X 10* 


1 


5.96 X 10^ 


1.6-2.3 


10 


R136 

It:'; i „i 


2-3 


10^ 


1.1 


8.97 X 10^ 


2.2-2.6 

4C4-^l4--, -,4- .,1 


11-13 



^Fieer et al. (1999) . ^Stolte et al. (2002) . ^Stolte e t al. (2005). *Stolte et al. (2006) . ' 'Kim et al. (2006) 
' ^Espinoza et al. (2009) . Niirnberger fc Petr-Gotz cns (200^, ^Sung fc Bessell (2004) . ^Haravama et al. (2008) 
l °Brandner et al. (2008) . "Andersen et al. (2009) . ^ ^Brandl et al. (1996), ^ ^Massev fc Hunter (1998) 



Galactic GCs 
UCD sim(SFE l.O.HE 1.00) 
UCD Sim (SFE 0.4,HE 1.00) 
UCD Sim (SFE 0.4.HE 0.03) 
-0.94»log||,(p„/10'Mopc-')+2.14 
NGC 3603 
Arches 
Wdl 
R136 



canonical IMF slope 



0.01 




3.5 



Galactic GCs 
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♦ 
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0.01 0.1 1 10 100 

Cloud density [10^ MqPc"^] 



1000 



Figure 2. The high-mass IMF slope, ag, for GCs (open circles 
with errorbars) decreases with increasing pre-GC cloud-core mass, 
Mci (stars-f gas). The solid line is a least-squares fit to the GC 
data and the corresponding equation is indicated in the legend. 
The IMF becomes canonical (03 = 2.3, horizontal dotted line) be- 
low 6.8 X IO^Mq (vertical dotted line). The range of quoted MF 
slope values found in the literature for the massive, young clus- 
ters NGC 3603, Arches, Wdl and R136 are indicated as the open 
pentagon, diamond, up- and downward triangle, respectively, with 
errorbars, being derived from the range of published 03 measure- 
ments. The overall trend is that more massive GCs form relatively 
more massive-stars (flatter IMFs). The filled symbols correspond 
to the initial conditions of UCD simulations that lead, after pri- 
marily gas expulsion and stellar mass loss driven evolution, to 
objects that resemble the properties of observed UCDs today. 
Different symbols correspond to different input parameters in the 
UCD m odels (star formation efficien cy SFE and heating efficiency 
HE, see iDabrinehausen et ai]|201Cll) . The UCD data form a sepa- 
rate group that runs roughly parallel to the GC data. Survival of 
GCs with 03 < 1 (grey-shading) is questionable (Sec.O. 



"^(^)^| 2.3, otherwise, ^^^^ 

A is either logio(A'/ci/lO^M0) (Fig. [U, 
logio(Pci/lO«M0) pc-3 (Fig. m or [Fe/H] (Fig. O 
with the corresponding coefficients, px and qx. The param- 
eter Aiini is the limiting value above or below which the 
IMF is top-heavy. The parameters are shown in Tab. [S] 



Figure 3. The high-mass IMF index, ag, versus the cloud-core 
density, p^l (stars+gas). Symbols and lines as in Fig. [2] The gap 
between the GCs and UCDs seen in Fig. |2]vanishes and the IMF 
indices, 03, now form a single sequence as a function of p^i: The 
IMF index at the high-mass end decreases (the IMF becomes 
flatter) with increasing density of the objects' progenitor cloud- 
core. It is canonical (03 = 2.3) below 9.5 X 10* Mq pc~^. 



Table 4. IMF indices, in equation ^ in dependence of the 
mctallicity, [Fe/H]. The indices «i/2 calculated from equa- 
tion I II2I I. and «3 is obtained using equation JTTJ with the param- 
eters in Tab. [3] For [Fe/H]= 0.0 the IMF is canonical everywhere. 
The IMF metallicity- variation is depicted in Fig. [5] 



[Fe/H] 


01 


«2 


as 


-2.0 


0.30 


1.30 


1.31 


-1.5 


0.55 


1.55 


1.64 


-1.0 


0.80 


1.80 


1.97 


0.0 


1.30 


2.30 


2.30 


+0.5 


1.55 


2.55 


2.30 


+ 1.0 


1.80 


2.80 


2.30 



with the appropriate inequality sign (< or >) for Aum being 
reported. 

The variation of the stellar IMF based on the metallic- 
ity is summarized in Fig. [S] Above a stellar mass of IMq 
the IMF flattens with decreasing metallicity (Fig. |4l equa- 
tion . By considering young and intermediate-age open 
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>&l Galactic GCs 

— a,=0.66"'[Fc/H]+2.63 

fNGC 3603 
Arches 
Wdl 
\r R136 



canonical IMF slope 




-2 -1.5 -1 -0.5 

global cluster metallicity [Fe/H] 

Figure 4. The high-mass IMF index, 03, as a function of the 
global cluster metallicity, [Fe/H]. Symbols and lines as in Figs. [2] 
and [3] NGC 3603, Arches and Wdl have about solar metallic- 
ity. For R136 [Fe/H]ai —0.48 is adopted, a typical value for the 
LMC. The data suggest the IMF to become less top-heavy with in- 
creasing cluster metallicity, having the canonical value (03 = 2.3) 
above [Fe/H] = -0.5. 



S' 
0) 
Q. 
O 



a3=-0.4072x'+1.9383 — 

NGC 3603 
Arches 

Wdl A 

R136 ^ 



canonical slope 




-0.1 405*[Fe/H]+logi o(Pc|/1 0^ Mq pc"^) 



Figure 6. Edge-on view of the fundamental plane, p, (equa- 
tion [I4l solid regression line) for the IMF index, Q!3, cluster den- 
sity, logj^QPcli a-nd metallicity, [Fe/H]. The abscissa is the new 
axis, x' , in which the scatter around the solid regression line 
through the data projected onto the a^—x' plane is minimal (equa- 
tion [I3ll. 




1 10 100 

m[Mo] 

Figure 5. Suggested shape of the stellar IMF for different metal- 
licities, [Fe/H] (not taking into account the density dependence 
of the IMF). The IMFs are scaled such that their values agree 
at m = IMq. Above IMq the IMF slope is determined by the 
present wor k (Fig. equa tion II 11 1. Below IMq the parametri- 
sation is by iKroupal 1120011 . equation I12I I , whose results suggest 
tentative evidence that more metal-rich environments produce 
relatively more low-mass stars. Note that only the metallicity de- 
pendence is shown, but not the dependence on mass (Fig. [2} or 
density (Fig. |3]l. 



clusters iKroupal (|200ll ) presented evidence that the IMF 
might also be dependent on [Fe/H] below IMq in the sense 
that a flatter low-mass IMF occurs for lower metallicities. 
His parametrisation for low-mass IMF variation is adopted 
here, 

Qi/2 = Qi/2,c + AQ[Fe/H] , (12) 

where Aa « 0.5 and ai/2,c are the respective slopes of the 
canonical IMF (Sec. The variation described by equa- 
tion [12] is incorporated in Fig. [5] Note that this parametri- 
sation has been suggested for metallicities > —0.5. For the 
lower metallicities needed here the calculated low-mass IMF 
slopes are thus extrapolated values. Furthermore this varia- 
tion is not incorporated in the lMarks fc Kroupal l|2010l ) mod- 
els and cannot be used to explain the observed low-mass star 
depletion in the observed GCs: Using equation 112! the IMF 
becomes steeper with increasing metallicity, contrary to ob- 
servations of the PDMFs (Sec.[T|. The resulting IMF indices, 
Qi, are reported in Tab. ID 



4 A FUNDAMENTAL PLANE FOR 03, logpci 
AND [FE/H] 

In this section a formula, 03 (log p^i, [Fe/H]), is derived in or- 
der to describe the dependence of as on density and metallic- 
ity for GCs simultaneously. The idea is to find a fundamental 
plane, p, in (0:3, log pci, [Fe/H] )-space which minimizes the 
scatter of a3-values around p. 

The method used to fi nd p follows the same basic princi- 
ple as lKroupa et al.l l|2010l ) used to determine the orientation 
of the disc of Milky Way satellites. The coordinate system 
{as, log pel, [Fe/H]) is rotated stepwise by At? — 1° around 
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the aa-axis. In each step the original data are projected onto 
the new a3-a;'(i9) plane, where 

x'{i)) = cosi9x [Fe/H] + sini9xlogio(pci/10'' Mq pc"^),(13) 

and a linear regression in this plane is performed. The ro- 
tation angle i9 in which the root mean squared-deviation of 
the projected data from the regression line along the aa-axis 
is minimal is the angle with the least scatter. For the data 
in Tab.[T]this is the case for — 98°, i.e. x' lies close to the 
original as-pd plane. 

Fig.^depicts x' (actually a;7sin(98°)) versus the IMF 
index Q3. The regression line through the data in this rep- 
resentation is the fundamental plane, p, seen edge-on (i.e. 
the projection of p onto the a^-x' plane). The fundamental 
plane equation thus reads 



aaix') 

-0.4072 xx' + 1.9383 x' ^ 0.87 
2.3 otherwise 



(14) 



and making use of equation (I13|l for i9 — 98° and x' ^ 0.87, 

Q3(log/9ci, [Fe/H]) = 0.0572 x [Fe/H] 

-0.4072 X logio(pci/10'' Mq pc"^) + 1.9383 . (15) 

As evident from its larger coefficient in equation (|15p . den- 
sity has a stronger influence on Q3 than metallicity. This is 
expected given the significantly lower scatter in Fig. |3] in 
comparison with Fig. [l] The fundamental plane has only a 
slightly smaller scatter in 03 than the Q3-logpci plane. The 
summed squared-residuals are 0.56 and 0.57, respectively. 

The stronger dependence on logpci can be understood 
as follows: A dense cluster forming from gas with solar 
metallicity will still produce a top-heavy IMF, although 
the [Fe/H] dependence alone (Fig.[4| suggests the canonical 
slope. Vice versa, a more extended cluster with low- [Fe/H] 
will not be as top-heavy as expected from Fig. U 

The procedure used here is similar, but not equivalent, 
to a principal component analysis. 



5 SURVIVAL OF GCS WITH TOP-HEAVY 

INITIAL MASS FUNCTIONS AND CAVEATS 

The combined information from this work, from the inde- 
pendent analysis of the top-heaviness in UCDs and the (rare 
and weak) observational evidence of top-heavy IMFs in star- 
burst clusters, which are remarkably consistent with each 
other, suggests that some presently old dense stellar sys- 
tems might have formed stars which were selected from a 
top-heavy IMF. The flatness of the high-mass end of the 
IMF reflects the metal-content and the density of the envi- 
ronment in which the respective objects were born. 

The question may be raised whether embedded clusters 
with extremely top-heavy IMFs (in particular 03 < 1) would 
survive stellar evolution and, if they do, whether they are 
long-lived, because a signiflcant amount of mass would be 
lost from massive evolving stars, triggering cluster expan- 
sion. Actually, for Q3 < 1 (grey-shaded area in Figs. [21 El HI 
and [6]) almost the total mass of the cluster is contained in 
stars with m > IMq and ~ 90 per cent of the total mass will 
be lost due to stellar evolution. In the best case this mass 
is lost adiabatically, i.e. slowly with respect to the crossing- 
time of the cluster. The cluster then expands by a factor of 



10 (equation [6]) through stellar evolution alone. The expan- 
sion due to stellar evolution becomes st ronger if the cluster is 
mass-segregated (Vespcrini e t al.ll2009l '). as the cluster mod- 
els in Marks fc Kr oupa (201(j) are. In Tab. [1] the number- 
and mass-fractions of stars more massive than IMq, 



and 



r 

JO. 



08Mq 



£^cijn)dm, 



Jo.o 



(16) 



(17) 



respectively, as wel l as the total mass lost d ue to stellar 
evolution f following iDabringhausen et al.ll2010l ) for the here 
estimated a3-values are reported. 

Liighausen et al. (in preparation) investigate, by means 
of A''-body computations, the survival of star clusters with 
a top-heavy IMF having stellar masses between 5000 and 
25000 Mq and half-mass radii between ~ 0.4 and 1 pc. 
They flnd that their models, having average stellar densi- 
ties within the half-mass radius up to « 2 x 10"* Mq pc"'', 
with an IMF index 03 1.0 immediately dissolve through 
strong stellar evolution driven mass loss and an acceler- 
ated dynamics (for a flxed cluster mass the relaxation-time 
becomes shorter with increasing top-heaviness due to the 
smaller number of stars). Slightly less top-heavy models are 
not very long-lived. The birth stellar densities for the Galac- 
tic GCs (Tab. [TJ are however two to four orders of magni- 
tudes larger than the density-range covered by Liighausen 
et al. (in preparation). In particular survival is likely to be- 
come possible when the crossing time, tcr, becomes shorter 
than the stellar evolution time-scale m/m for massive stars, 
where m is the mass-loss rate of a star of mass m. There- 
fore survival of such dense objects with extremely top-heavy 
IMFs remains to be checked numerically. 

The constraints on the top-heaviness were derived as- 
suming the energy for gas removal is provided by stellar 
winds and radiation only and all energy radiated by the 
stars is completely ab sorbed by the ga s. This is likely to be 
an oversimplification IjPale et al.ll2005f ). If the initially most 
massive clusters are able to retain the residual-gas in their 
deep potentials until supernovae occur, their energy input 
might be sufficient to remove the gas in ~ 10^ yr. Addi- 
tionaly, the metal content might influence the gas throw- 
out through radiation-f-^gas/dust c oupling leading to differ - 
ent energy absorbtion efficiencies (|Marks fc Kroup"all201Gl ). 
Both effects need to be considered in future work. 

The parametrisation (equations [Tl] and I14p should 
therefore not be seen as an established dependence on the 
cluster initial parameters. Instead it is suggestive of what 
IMF variation the present observational data suggest. 



6 CONCLUSIONS 

This work has shown that in oder to remove the residual- 
gas within ~ 10^ yr after star formation from p rogeni- 
tors of some Galactic GCs (|Marks fc Kroupall2010l ') start- 
ing with a universal low-mass IMF, top-heavy IMFs need 
to be invoked. The IMF index/slope Q3 in equation ^ 
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takes values between 1 (very top-heavy) and 2.3 (the canoni- 
cal Salpeter/Massey- value). Although long-term survival for 
some GCs with strongly top-heavy IMFs might be at threat 
through stellar evolution alone, self-consistent modelling of 
very dense GCs with a top-heavy IMF is needed to inves- 
tigate their survival chances. The index 03 anti-correlates 
roughly log-linearly with the pre-gas expulsion cloud-core 
mass and density (Figs. [2] and [3]) and correlates linearly with 
the GC metallicity (Fig.gl. 

Combining the density and metallicity dependence of 
Q3, a fundamental plane has been found which describes 
the variation of the IMF with both parameters simultane- 
ously. The fundamental plane suggests that as varies more 
strongly with density than with metallicity. 

Large mass-to-light ratios in UCDs indicate this 
class of o bjects to have formed wit h a top-heavy IMF 
as well ( Dabringhausen et al.l I2OO8I ). The data from 
iDabringhausen et al.l (|20ld f on probable initial conditions 
of observed UCDs form a separate group of objects in the 
Q3— cloud-core mass diagram but run roughly parallelly to 
the GC data (Fig. [2]). The disappearance of this gap when 
Q3 is depicted against the cloud-core density (Fig. ^ might 
suggest th a,t UCDs form from mergers of GCs in cluster 
complexes (iFellhauer fc Kroupall2005l : lAssmann" et al.ll201ll : 
iBriins et al.ll201ll '). which moves merged GCs of the same 
top-heaviness to higher UCD masses in Fig. [2] 

The decreasing trend of 03 with increasing density is 
qualitatively understood since denser systems are expected 
to form more massive stars with respect to the canonical 
IMF, if a channel for massive-star formation is accretion- 
driven coagulati on of proto-stellar cor es due to a higher colli- 
sion probability (jMurrav fc Linll996h . And in systems of low 
metallicity the Jeans mass is expected to be larger favouring 
the formation of more high-mass stars (|Adams fc Fatuzzq 
I199& : Larson 1998). 

Are these dependencies compatible with observations 
of young clusters in the Local Group? For the MW and for 
the LMC the answer is positive (Figs. [2] |3] HI . Although 
the known Local Group starburst clusters lie just in the 
density-regime where the IMF starts to become top-heavy 
{pci ^ 9.5 X W'^Mq pc~^) the observations are consistent 
with the results. For the densities of these clusters, stellar 
dynamical biases would probably hid e any true present top- 
heavy IMF due to m ass segregation (Por tegies Zwart et al.l 
l2007f i. lMassevl l|2003l ) shows homogeneously estimated mass 
function slopes for OB associations in the MW and LMC 
finding all of them to be consistent with the Salpeter slope 
(q = 2.35). For the LMC clusters ([Fe/H]?5i -0.4) and solar 
metallicity MW clusters this is consistent with the results 
in this work (Fig. |4|, given the star-forming densities in the 
Local Group are < 10*^ Mq pc~^. 

A similar, but slightly different parametrisation of 0:3 
with the stellar cluster mass, M^ci , (Tab. [3]) has been 
made use of in IWeidner et al.l (|201ll 'l to include the ev- 
idence for a systematically varying IMF into the frame- 
work of the theory of the integrated ga lactic IMF (IGIMF, 
i.e. t h e stellar IMF of wh ole galaxies, iKroupa fc Weidneij 
l2003l : IWeidner fc Kroupal I2OO5). Under conditions when 
the galaxy- wide star formation rate SFR> 27 Mq yr~^ 
star burst clusters with masses Meci > 10^ Mq form 
ijWeidner et al.1 l2004h . The IGIMF, being the sum of all 
IMFs in all young star clusters, then also becomes top- 



heavjQ. Their models lead to reasonable agreement with 
cosmological observations. The results obtained here are 
cxplicitcly incorportated into the IGIMF framework in 
JKrmipa et al. (2 012| ), where, for example, it is shown that 
the IGIMF of the Galactic bulge agrees with the chemical 
abundance constraints. In which way the chemical evolution 
of galaxies over cosmological time would be affected by a 
top-heavy IGIMF is subject to further studies. 

To conclude, varying the low-mass IMF in order to ex- 
plain the observed PDMFs (and, in turn, to avoid top-heavy 
IMFs) in GCs does not appear to be a feasible solution. GCs 
formed with the canonical IMF over all stellar masses and 
evolving secularly over a Hubble time cannot lead to the 
observed concentration-PDMF correlation, unless they are 
born mass-segregated and filling their tidal radii. But this 
begs the question how such distended mass-segregated clus- 
ters can form in the first place. 

Instead, initially concentrated, dense < 1 pc large 
mass-segregated massive progenitors of present-day GCs 
uncover evidence for a systematically varying high-mass 
IMF with density and metallicity, in remarkably good 
agreement with the indepen dently obtained evidence fo r 
top-heavy IMFs in UCDs by IDabringhausen et al] (|2010t ). 
The suggested trend is qualitatively in agreement with 
the expectation that higher density and lower metallicity 
environments should form more top-heavy IMFs (see 
Sec. [T]). Observations are consistent with the here derived 
trend but mass function measurements for initially denser 
clusters than NGC 3603 and the Arches cluster are likely 
needed to conclusively test the here suggested behaviour 
for the high-mass IMF. 
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